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We examine quenched chiral logarithms in lattice QCD with overlap Dirac quarks. From our data of m^, 
we determine the coefficient of quenched chiral logarithm 5 = 0.203(14), 0.176(17), 0.193(17) and 0.200(13) for 
lattices of sizes 8^ x 24, 10^ x 24, 12^ x 24 and 16^ x 32 respectively. Also, for the first three lattice sizes, we 
measure the index susceptibility of the overlap Dirac operator, and use the exact relation between the index 
susceptibility and the rj' mass in quenched chiral perturbation theory to obtain an independent determination of 
d = 0.198(27), 0.173(24), 0.169(22), which are in good agreement with those determined from m^. 



1. Introduction 

In quenched chiral perturbation theory (qxPT) 
jl], the pion mass to one-loop order reads as 

ml = Cm - S[ln{Cmg/Al) + 1]} + Bm^ (1) 

where denotes the bare {u and d) quark mass, 
is the chiral cutoff which can be taken to be 
2V27r/7r ( fir — 132 MeV ), C and B are parame- 
ters, and 6 is the coefficient of the quenched chiral 
logarithm. 

Theoretically, S can be estimated to be p| 



S = 



(2) 



where m^/ denotes the 77' mass in qxPT, and Nf 
is the number of light quark flavors. For = 132 
MeV, iV/ = 3, and m^, = 
853 MeV, (|) gives 



2m\ = 



,5-0.176 



(3) 



Evidently, if one can extract 5 from the data of 
771^, then 771,,' in qxPT can be determined by (||). 

Besides from the data of 777,^, one can also ob- 
tain 5 via d^) by extracting 777,,,/ from the propa- 
gator of the disconnected hairpin diagram. How- 
ever, to compute the propagator of the hairpin is 
very tedious. 
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Fortunately, with the realization of exact chiral 
symmetry on the lattice |^,^ , the quark propaga- 
tor coupling to 77' is {Dc + mq)~^ ||], thus only 
the zero modes of Dc can contribute to the hair- 
pin diagram, regardless of the bare quark mass 
mq . Therefore one can derive an exact relation ||] 
between the 77' mass in qxPT and the index sus- 
ceptibility of any Ginsparg- Wilson lattice Dirac 
operator, without computing the hairpin diagram 
at all. Explicitly, this exact relation reads as 



{rrin 



ANf ((77+-77_)2) 



(4) 



where Ng is the total number of sites, and x = 
((77+ — n_y)/Ns is the index susceptibility of 
any Ginsparg- Wilson lattice Dirac operator in the 
quenched approximation. Then (^) and to- 
gether gives 



1 



((7 



2^2(/,a)4 



Ng 



(5) 



A salient feature of (^ is that S can be determined 
at finite lattice spacing a, by measuring the index 
(susceptibility) of the overlap Dirac operator, and 
with /ttA extracted from the pion propagator. 

Now it is clear that, in order to confirm the 
presence of quenched chiral logarithm in lattice 
QCD, one needs to check whether the coefficient 
S obtained by fitting to the data of 777^ , agrees 
with that (|) fr om the index susceptibility. This 
is a requirement for the consistency of the theory, 
since the quenched chiral logarithm in ml (|l|) is 
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Table 1 



Determination of the coefRcicnt of quenched chiral logarithm 6 by fitting (Q) to data of m^. The ij' mass 
in qxPT is computed from (^). 









(MeV) 














xT 


a(fm) 


min. 


max. 




a/ A™" 


Ca 


B 


S 


m^/ (MeV) 


8^ X 24 


0.147(1) 


418 


1047 


2.5 


0.78 


1.193(18) 


1.152(56) 


0.203(14) 


916(32) 


10^ X 24 


0.152(1) 


375 


834 


2.9 


0.64 


1.017(25) 


1.222(71) 


0.176(17) 


853(42) 


12^ X 24 


0.147(1) 


395 


870 


3.5 


0.65 


1.051(24) 


1.254(73) 


0.193(17) 


893(40) 


16^ X 32 


0.095(1) 


439 


1101 


3.4 


0.53 


1.117(15) 


2.146(102) 


0.200(13) 


909(31) 



due to the rj' loop coupling to the pion propagator 
through the mass term (in the chiral lagrangian), 
thus S (to,,') must be the same in both cases. We 
regard this consistency requirement as a basic cri- 
terion for lattice QCD (with any fermion scheme) 
to realize QCD chiral dynamics in continuum. 

Recently, we have determined S = 0.203(14), 
from the data of to^, as well a.s 6 — 0.197(27) 
from the index susceptibility, for the 8'^ x 24 lat- 
tice 1^. Their excellent agreement suggests that 
lattice QCD with overlap Dirac quarks indeed re- 
alizes the QCD chiral dynamics in the continuum. 

In this paper, we extend our studies to other 
lattices of larger volumes and/or smaller lattice 
spacings. Our computations are performed with 
a Linux PC cluster of 30 nodes at the speed ~ 30 
Gflops §. 

2. Results 

With the Wilson SU{3) gauge action and 
Creutz-Cabbibo-Marinari heat bath algorithm, 
we generate ensembles of gauge configurations as 
foUowsR. 



L '^ xT 


/3 


# configs. 


8^ X 24 


5.8 


100 


10^ X 24 


5.8 


100 


12^ X 24 


5.8 


100 


163 X 32 


6.0 


56 



For each configuration, quark propagators are 
computed for 10-12 bare quark masses. Then the 
pion propagator and its time correlation function 
G{t) are obtained, and {G{t)) is fitted by the 



usual formula 



l^-m^at ^ g-m,a(T-t)l 



2TO^a 

to extract the pion mass m^a and the pion decay 
constant 



f^a = 2mqa- 



^The 56 configurations on t he 16'^ x 32 1 attire are rp tripvRd 
from the gauge connection ( ittp://qcd. nersc.gov/ ). 



The data of /jra can be fitted by straight line. 
Thus taking /^a at rUqa = equal to 132 MeV 
times the lattice spacing a, then a can be deter- 
mined. Fixing A-^a = 2\/27r/^a, we fit ^ to the 
data of (m^a)^, and obtain Ca, 6 and B. Then 
to^ is computed from (^). Our results are sum- 
marized in Table [|. Evidently, all 6 values are 
in good agreement with the theoretical estimate 
6 ~ 0.176 from q^PT. Furthermore, the r]' masses 
are also in good agreement with the theoretical 
estimate to,,' ~ 853 MeV. 

Even though the quenched chiral logarithm 
may not be easily detected in the graph of 
(mTrfl)^ vs. TOqfl, it can be unveiled by plot- 
ting (TO7ra)^/(m^a) vs. m^a, as shown in Fig. 
|l|. Further, plotting (TO^a)^/(TOqa) — B{mqa) vs. 
log(TOqa), the presence of quenched chiral loga- 
rithm is evident, as shown in Fig. |^. 

Note that, in Fig. 0, for the 16^ x 32 lattice at 
l3 = 6.0, the minimum of (jn^a)'^ /{niqa) occurs at 
rriqa = Ca6/B ~ 0.1, which corresponds to rriq ~ 
200 MeV for a ~ 0.095 fm (Table 0). Similarly, 
for 8^ X 24, 10^ x 24, and 12^ x 24 lattices at = 
5.8, the minima of {rriT^af' / {mqa) occur at m^a ~ 
0.21, 0.15 and 0.16 (i.e., to, ~ 280, 190, and 217 
MeV) respectively. Now, except for the 8^ x 24 
lattice which has the smallest volume, the minima 
of {m.^af' / {rnqo) for the other three lattice sizes 
almost occur at the same to„ ~ 200 ± 20 MeV. 
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Table 2 

Determination of S (|^) from index susceptibility. 
The 77' mass in qxPT is computed from (W). 



(m^ a)^ / (mq a) - B ( rriq a ) 
1.6 





X T 


x/10-' 


6 


m,,/(MeV) 


8^ 


X 24 


3.67(50) 


0.198(27) 


904(61) 


10^ 


X 24 


3.65(50) 


0.173(24) 


817(58) 


123 


X 24 


3.13(41) 


0.169(22) 


835(55) 



This provides rather substantial evidence that the 
quenched chiral logarithm parameters extracted 
from our data of is indeed genuine, not due 
to any finite size effects or discretization errors. 



( a )'/ ( mp a ) 
1.55 T 



16 X 32 
Wilson, p = 6.0 
56 configs. 




.J 2Q 20D2/7/8 twqcd/l1Bt32/b6O/m013/g56/mpi2omq 
0.00 0.05 0.10 0.15 0.20 0.25 

mq a 

Figure 1. {rriT^a)'^ /(mqa) versus the bare quark 
mass m„a. 



Next we measure the index susceptibility of 
the overlap Dirac operator by the spectral flow 
method |3|, and then determine 5 via The 77' 
mass in qxPT is computed from (Q). Our results 
are summarized in Table including our earlier 
results for the 8'^ x 24 lattice. Evidently, each 5 is 
in good agreement with the corresponding S ex- 
tracted from toJ (Table Q), as well as with the 
theoretic estimate 0.176. 
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16^ X 32 
Wilson, p =6.0 
\ 56 configs. 




V 8 
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= 0.200(13) 
= 2.146(102) 


twqcd/l ■ 


6t32/b60/m013/g56/chilog 


0.1 




log( mq a ) 





Figure 2. The extraction of the quenched chiral 
logarithm by plotting {m.^af' / [mqa) — B{raqa) 
versus log(m^a). 



3. Conclusion 

Our results in this paper provide strong evi- 
dences that quenched QCD with overlap Dirac 
quarks realizes quenched QCD chiral dynamics, 
as depicted by quenched chiral perturbation the- 
ory 
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